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We calculate single atom heating rates in a far detuned optical lattice, in connection with recent 
experiments. We first derive a master equation, including a realistic atomic internal structure and 
a quantum treatment of the atomic motion in the lattice. The experimental feature that optical 
lattices are obtained by superimposing laser standing waves of different frequencies is also included, 
which leads to a micromotional correction to the light shift that we evaluate. We then calculate, and 
compare to experimental results, two heating rates, the "total" heating rate (corresponding to the 
increase of the total mechanical energy of the atom in the lattice), and the ground bande heating 
rate (corresponding to the increase of energy within the ground energy band of the lattice). 
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I. INTRODUCTION 

In the field of atomic quantum gases, optical lattices 
have become a versatile tool to trap atoms in an almost 
non-dissipative way. This allows one to simulate many- 
body Hamiltonians originally formulated in condensed 
matter physics, such as the Hubbard Hamiltonian for 
bosons and fermions (see [l[ for a recent review). This 
also opens promising implementations of quantum logi- 
cal operations 0. For all these applications, decoherence 
due to residual spontaneous emission of the atoms in the 
optical lattice has to be kept as low as possible. It was 
realized recently, on an experimental implementation of 
the Bose Hubbard model, that a noticeable heating rate 
of the atoms takes place and has to be included to obtain 
fair agreement with the theory 

In this article, we perform a theoretical analysis of the 
increase rate of the atomic mechanical energy in a far- 
detuned optical lattice. Contrarily to earlier laser cooling 
studies, relying on semiclassical approximations 0-0], or 
restricting to the Lamb-Dicke limit [Io|, U3 j or consider- 
ing reduced dimensionalities and simplified level schemes 
Ell EH j we include the relevant fine and hyperfine 
atomic structure and fully take into account the quan- 
tum motion and tunneling of the atoms in the three- 
dimensional lattice. We however restrict for simplicity to 
the single atom problem: We thus miss the photoasso- 
ciation processes, which mainly produce losses of atoms 
[Til ], and the extra heating due to multiple scattering of 
fluorescence photons in the atomic gas. Fortunately, as 
predicted in [15|, [l|| and observed experimentally in [l7| , 
this extra heating is small in the so-called festina lente 
limit, where the fluorescence rate is much smaller than 
the oscillation frequency of an atom in the lattice, a re- 
alistic regime for far-detuned optical lattices. 

The article is organized as follows. We present our 
model for the atomic structure and for the laser field 
producing the lattice, and we derive a master equation 
for the ground state atomic density operator in section 
HIl We define and calculate two types of heating rates in 
sections IIIII and IIV1 In section IIIII we show analytically 




FIG. 1: The atomic level scheme considered in this work, 
typical for alkali atoms. 

that the total heating rate, defined as the time deriva- 
tive of the total mean mechanical energy of the atom, 
remarkably is independent of the atomic quantum state 
and insensitive to the sign, blue or red, of the laser de- 
tuning with respect to the atomic resonance. In section 
HVl we calculate the ground band heating rate, that is 
the increase rate of the energy within the ground Bloch 
band of the lattice, for an atom initially in the ground 
Bloch state; we show that it widely depends on the blue 
or red nature of the lattice in the tight-binding limit. We 
conclude in section fVl 



II. MODEL AND GROUND STATE MASTER 
EQUATION 

A. Atomic structure 

We consider in this article alkali atoms, with an in- 
ternal level structure shown schematically in Fig. [TJ We 
consider only the ground s state and the first p excited 
states, and neglect the other excited states in our anal- 
ysis. We denote by g the ground state manifold, and by 
ej, with J — 1/2 or 3/2 labelling the electronic angular 
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momentum, the two excited fine multiplets (leading to 
the so-called Z?i and D 2 lines) separated in energy by an 
amount fkuj from the ground state. The fine structure 
atomic levels are further split by the hyperfine interac- 
tions -ffhf, giving rise to hyperfine multiplets with total 
angular momentum F = I ±1/2 in the ground s state (de- 
noted by 31 and ff2 ) and F' = 7±l/2, F 1 = J±3/2,/±l/2 
in the e\/2 and 63/2 manifolds, respectively. For order of 

magnitude estimates to come, we will note Ay C typical 
values for the hyperfine splittings in the ground or ex- 
cited states, respectively. The atomic density operator 
tr has thus many internal atomic matrix elements, that 
may be collected in ground state elements, excited state 
elements and optical coherences, forming the submatrices 

Ogg = Pg(TPg, O e j e ; , = P e ,<jP ej ,, '&YL& <J e j g = P e , (7 Pg . 

Here P g projects onto the ground state manifold (includ- 
ing the hyperfine splitting in the two ground states g\ and 
32); and P ej projects onto the excited state fine multiplct 
ej- 

B. Laser configuration 

In the experiment of Q, a three-dimensional cubic 
optical lattice was created by superimposing standing 
waves along x, y, z, with orthogonal linear polariza- 
tions and with widely different frequencies. The rapidly 
oscillating interference terms in the laser intensity be- 
tween standing waves along orthogonal directions aver- 
age to zero and the resulting light shift potential is es- 
sentially scalar. We include this multichromatic feature 
in our model: Writing the total laser field as a sum of 
positive and negative frequency components, E(r,t) = 
E(+)(r,t)e~^ I - t + E^)(r,t)e i " I - t , with ui L being the car- 
rier frequency, the amplitudes E^ (r, t) are taken as 

E(+)(r,i)= £ e^(r)e-" A ^ (1) 

fj.=x,y,z 

and E(")(r,t) = [E(+)(r,t)]*, e M being the unit vector 
along direction [i £ {x, y, z}. The modulation frequencies 
A^ associated with each axis n = x,y,z are assumed to 
be incommensurate and much smaller than the carrier 
frequency ojl- We define the detunings from the excited 
states as 

Sj =ul— (2) 

where ujj is the central resonance frequency of the atom 
on the g — >• ej transition, see Fig. [1] For simplicity, we 
assume that these detunings are much smaller than the 
atomic resonance frequency, so that the atom-laser cou- 
pling can be taken in the rotating wave approximation 
(KWA) . This is consistent with our approximation where 
only the dominant s — p transition is considered. How- 
ever, the detunings are assumed larger than any other 
frequency scale in the problem, including the residual 
modulation frequencies 

A mod « \S\. (3) 



Here and in what follows, when estimating orders of mag- 
nitude, we will write for simplicity 5 as a typical order of 
magnitude for the detunings 5j and A mo d for the modu- 
lation frequencies A^. 

In the rotating frame, the atom-laser coupling in the 
rotating wave approximation reads 

V AL (t) = -D<+) • #)(r, t) + h.c. (4) 

where D( + ) is the raising part of the atomic dipole op- 
erator. The resulting time-averaged light shift potential 
is almost scalar, since, for the considered atomic struc- 
ture, the following ground state "polarizability" operator 
is scalar, 

{T>-B L ) gej {T>-e L ) ejg = d 2 jP g (5) 

where is any unit vector with real components, and 
(X) ejg = P ej XP g is the restriction of any operator 
X between the manifolds ej and g. Eq. ([5]) may be 
checked in the fine structure basis from the expression 
of the Clebsch-Gordan coefficients, where it appears as 
a well known property of 1/2 — > 1/2 and 1/2 — > 3/2 
transitions. From the orbital rotational invariance of 
(D) eff • (D) fle = d 2 P e , one also deduces the relation 

( D )e J9 ' ( D ) 9 e. ; , = d'Pe., Sjj, (6) 

where d is the so-called reduced atomic dipole moment 
for the s — p transition. This finally leads to 

d 2 2d 2 

d l/2 = Y ^3/2 = — ■ ( 7 ) 

C. Equations of motion 

The starting point of our treatment is the fully quan- 
tum optical Bloch equation for the atomic density opera- 
tor a [18| , including the hyperfine structure and a quan- 
tum treatment of the external atomic variables (center of 
mass position). We assume the following, typical hierar- 
chy between the different energy scales in the problem: 

E Iec , E kin , U max < HT < \hS\. (8) 

Here E Tec is the atomic recoil energy, 

Erec = (9) 

2m 

k,L = i^l/c is the laser wavevector, T is the natural 
linewidth of the excited states, -Ekin is the atomic ki- 
netic energy, and C/ max ~ \hfl 2 /S\ is the typical depth of 
the optical lattice potential, with a laser induced Rabi 
frequency f2 loosely defined by Val ~ t&l- The con- 
dition t/ max <C h\S\ implies the weak saturation limit, 
o 2 

w <C 1. In these conditions, we can adiabatically elimi- 
nate the fast optical coherences and excited state matrix 
elements, which track the slowly-evolving ground state 
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variables. The large detuning regime Eqs. (|3l8|) consider- 
ably simplifies this elimination, allowing one to perform 
an expansion of a eg and a ee in powers of 1/6, up to order 
1/5 2 here. 

As detailed in the Appendix [XJ this treatment leads to 
an effective equation for the ground state density matrix 
o gg , which is still rather complicated due to the hypcrfinc 
Hamiltonian and to the explicit time-dependence intro- 
duced by the temporal modulation of the laser electric 
field at frequencies A M . The light shifts in particular 
have a non-scalar part, that also provides a Raman cou- 
pling between g\ and g2 with a Rabi frequency ~ fl 2 /S. 
The equations can be further simplified in the experi- 
mentally relevant case, where the modulation frequency 
is much faster than the atomic motion in the optical lat- 
tice, but much smaller than the ground state hypcrfinc 
splitting (so as to make the hyperfine Raman coupling 
widely non-resonant), 



fi 2 
S 



< A mod < A 9 



hp 



(10) 



The inequalities in Eq. (|10|) have two consequences. 

First, due to the smallness of the laser induced Ra- 
man couplings as compared to the ground state hyperfine 
splitting, ground state hyperfine coherences will depart 
from their initial zero value by a small amount, of order 
SI 2 /(SA^). Their adiabatic elimination (valid in the ab- 
sence of Raman resonances) leads to a small correction 
to the light shift of order 



SU] 



hf 



nn 4 



(ii) 



After this second adiabatic elimination, the unknown 
part of the density operator is now the diagonal part 
<jj?' as of <jgg, the diagonal part of any ground state oper- 



ator X being defined as 



X di ^ = P gi XP gi 



Pg 2 XPg 2 , 



(12) 



where P gi projects onto the ground state gi, i £ {1,2}. 
Of course, cr^ ag still contains Zeeman coherences within 
the g\ and gi manifolds. 

Second, due to the time-depcndcncc of the laser field 
at frequencies A M , er^ 6 can be decomposed as a slowly- 
evolving, zero frequency component cr^g 3 " 6 plus fast oscil- 
lating components at harmonics of the modulation fre- 
quencies A M . The latter correspond to an atomic micro- 
motion in the time-dependent optical lattice, similar to 
the dynamics in Paul- type ion traps (l9j . The fast com- 
ponents are typically much smaller than the slow one by 
a factor rj fl 2 / (SA moc \). In this regime, a perturbative 
technique [2(], Hl[ allows to construct a time-independent 
effective potential induced on <J gg as by this micromotion, 
which constitutes another small correction to the light 
shift, explicitly calculated in the Appendix [S] and of or- 
der of magnitude 



SUmic 



S 2 A 2 , 

mod 



(13) 



If A^j 3> E Tec A^ { , as is the case in Q, the hyperfine 
contribution in Eq. (jllj) dominates over the micromotion 
contribution in Eq. (|13j) . 

The hypcrfinc (jTTJ) and micromotion (fl~3"|) contributions 
appear at order S" 1 of the adiabatic elimination of a ee 
and a eg , and originate from a purely conservative (though 
time-dependent) light shift potential. At order 5~ 2 in 
the adiabatic elimination, one obtains non-conservative 
terms, proportional to the fluorescence rate r(fi/<5) 2 , and 
corrections to the light shift potential of order 



SU l 



/6* 



frA e ' 9 

nAA hf i 



KAr 



(14) 



All these terms may be directly time averaged, neglecting 
micromotion corrections at this order. 

This procedure results in the final master equation for 
the zero frequency hypcrfinc diagonal part of the ground 
state density operator, 



_^_-diag 

dt 99 



^2L + u + m ,^-± {Wtd %* } 



i 

2 



(15) 



The structure of this equation is familiar from earlier 
studies on laser cooling. The first commutator corre- 
sponds to a Hamiltonian evolution in the light shift po- 
tential U + SU. For our choice of polarizations and de- 
tunings, the basic scalar light shift potential U is scalar, 




Em* 



(16) 



The quantity SU, whose complete expression is given in 
the Appendix [2] includes all previously discussed cor- 
rections in Eq. (|11I13I14|) . As discussed in the caption 
of Table HI all these corrections are negligible to a good 
accuracy for typical experimental parameters. 

The part in Eq. (|15[) involving an anti-commutator 
with the operator W corresponds formally to an anti- 
hermitian Hamiltonian, and describes departure from the 
ground state upon absorption of a laser photon. This 
part is also scalar as expected, 



W(r)= E r 



d 2 j 



(frSj) 2 



P, 



£lMr)l 



(17) 



Finally, the last "feeding" term in Eq. (|T5| describes 
atoms returning to the ground state after an absorption- 
spontaneous emission cycle. This part involves a Liou- 
villian operator C acting on the excited state density op- 
erator, 

J,J' J e±n 

e -ife n.r / D (-) . £ *\ a / D (+) . A e tk n-v 

\ ) gej ' ' \ J ejig 

(18) 
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where r = d 2 / (3Treoh) is the usual spontaneous 
linewidth, and where the fine energy splitting in the mo- 
mentum fco of the spontaneously emitted photon was ne- 
glected. In ([15]) C acts on the zero frequency component 
of the excited state density operator expressed in terms 
of a aa as 



1 



diag 
fiJejg^n'J gg 



S.aTHD^e^;. (19) 



If one takes the trace over the internal atomic states of 
this expression, in the case J = J', one obtains from 
Eq. (0 the useful property 



Tri nt (<7 ejej ) 



-^^(r)Tr int (a^)^(r). 



III. TOTAL HEATING RATE 



From the master equation Eq. (|15p . we now calculate 
the rate of change dE to t /dt of the total mean mechanical 
energy of the atom: 



Etot = Tr[( 



2m 



U 



gg 



where averaging is done on both internal and external 
degrees of freedom. The first commutator in the right 
hand side of Eq. (fl~5"|) does not contribute to dEtot/dt. For 
the other two terms in the right hand side of Eq. (fl"5|) . 
we calculate exactly the contribution involving p 2 /2m 
and U in the energy, and we put a simple bound on the 
contribution involving SU. Since p 2 /2m and U are scalar, 
all internal traces may be evaluated, and we obtain 



d F 
dt 



I h 2 k 2 

j) 2 ) E (i^fi^wi 

+ 0[T(n/5) 2 5U}, (22) 

where (X) = Tr[X Vg 1 ^] is the expectation value of any 
ground state observable X. To obtain this result, we have 
taken the sum over e and the integral over n in Eq. (fT5|) , 
using e^on.r^g-tfcon.r = Hq _ p . nk n/m + h 2 kf/(2m), 
where Ho = p 2 /(2m) + U + SU; the term linear in n is 
odd and vanishes after integration over n. We also used 
Eq. @ and Eq. (gDJ). 

The last term in Eq. (|22[) , originating from SU , is neg- 
ligible as compared to first term, originating from the 
recoil due to spontaneous emission, provided that 



\SU\ < 



h 2 k 2 

2m 



(23) 



Using the estimates given in Eqs. (lllll3ll4[) we find that 



this is extremely well obeyed in [3j. Condition Eq. 
is supposed in what follows to be satisfied. 



Eqs.jn 
in [|~ 



Eq. (|22p can be further transformed, using the fact that 
Maxwell's equation imposes (A r + k\)£^(r) = where 
kh = lul/c and A r is the Laplacian operator. Then 



dt 



-E tt 



(h8j) 



2m 



|gradf M (r) 



(24) 



where the laser field current for the field amplitude £ M 
is iii = fr(£*grad£ M — c.c.)/(2im) and one also has 
|grad£ M | 2 = kHS^ 2 + A r \£ tJl \ 2 /2. Furthermore, if each 
(20) amplitude £ M is a simple laser standing wave, £^(r) = 
£o,fi cos(k M • r + ipfj,), and under the reasonable assump- 
tion that fco may be identified to the laser wavevector fc^, 
one finally obtains 



dt 



Etnt = 




(25) 



This may be rewritten in terms of the maximal fluores- 



cence rate in the lattice, that is the maximum of 



(21) W: 



21.2 



dt 



-E t , 



K 2 k 



2m 



rmax 
-i fluo ' 



(26) 



In this form, (|2"6"|) reproduces in the classical limit the 
position-independent value of the momentum diffusion of 
a two-level atom in a weak laser standing wave, obtained 
from laser cooling theory (29J. 

Eq. (f2~51) is one of our main results. It shows the coun- 
terintuitive fact that for a given laser intensity distribu- 
tion, the total atomic heating rate in a far-detuned op- 
tical lattice is independent of the atomic external state 
and of the sign of the atom-laser detuning, provided that 
each laser standing wave is linearly polarized and the 
atomic kinetic energy remains small as compared to HT. 
In particular, trapping the atoms at the nodes of a blue 
detuned optical lattice, although it reduces the atomic 
fluorescence rate, does not reduce the total heating rate 
dEtot/dt with respect to trapping at the antinodes of a 
red detuned optical lattice with the same absolute value 
of the detuning and laser intensity. 



IV. GROUND BAND HEATING RATE 

We now consider the increase rate of energy within 
the ground band of the optical lattice. For experiments 
with atomic gases, the physical motivation to do so is 
twofold: One can imagine observable quantities that 
depend mainly on the probability amplitudes to find 
the atoms in the Bloch states of the ground band, and 
one can also imagine that evaporative-type experimental 
techniques are developed to continuously eliminate the 
atoms from the excited energy bands. 
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FIG. 2: Sketch of the different heating processes in a lattice in 
the tight-binding limit: intraband processes (1) and interband 
processes (2). 



Having previously bounded the effect of the light shift 
correction 5U, we directly neglect it in the master equa- 
tion Eq. (|15j) . and we define the ground band mean en- 
ergy as 



Egb = 



Tr[P GB (£- + U)af^ 



gg 



Tr[P GB a 9 t S ] 



(27) 



where Pbg projects onto the ground energy band in the 
periodic potential U. The increase rate of Pgb, con- 
trarily to the total energy increase rate, depends on 
the atomic state. To make the analytical calculation 
tractable [3(| , we thus restrict to the initial increase rate 
^-Egb(£ = 0) for an initial atomic state only populat- 
ing the ground Bloch state, that is the state with zero 
Bloch vector q = in the ground band (of band indices 
0). This choice is motivated by the fact that the heating 
rate of [j| was measured for a Bose-Einstein condensate. 

We shall present numerical results in the general case, 
and then focus on the tight-binding limit, with a modula- 
tion depth of the optical potential U (r) much larger than 
the atomic recoil energy, where analytical results are ob- 
tained. It is then convenient to assume that the optical 
lattice potential has a minimum at the origin r = of 
the coordinates. We thus take £ M = 5o,m cos (^ ' r ) f° r a 
red detuned lattice (5j < VJ), and £ M = £q^ sin(k M • r) 
for a blue detuned lattice (Sj > VJ). 
General results: The initial atomic density operator is: 



<i aS (* = 0) = |O;q = O)<O;q = O| 



T int 



(28) 



Here cr™* is any ground internal atomic state (with no 
hyperfme coherences), |0;q) is the ground band eigen- 
state of Bloch vector q, normalized in an arbitrarily 
large quantization volume commensurate with the lat- 
tice spacing, and (for futur reference) Prj;q is the corre- 
sponding eigenenergy. Since 5gg as (i = 0) commutes with 
H = p 2 /(2m) + U(r) and with the projector Pgb, the 
time derivative of the ground band energy at t = has, 
from (|15p . the simple expression involving the feeding 



term only: 
dE G 



dt 



(f = 0) = Tr[P GB (Po - £o ;q =o)£[a ee (0)]. (29) 



Replacing the feeding operator by (|18j) and the excited 
state density operator by its expression (|19p . we obtain 
at time t = 0: 



dt 



= r 




\ J - -' / - fx 

x(£ ; o;fe L (e, l -n)-^o;o)|(0;fcL(e Al -n)|e- 



**-'£ M (r)|0;0)| 5 



(30) 



We have used the fact that, after absorption of a laser 
photon in the mode £ M and spontaneous emission of a 
photon of wavevector fe^n, where k was identified to 
kL, the initial center of mass atomic state of zero Bloch 
vector acquires a Bloch vector q = fc^e^ — k^n. The 
quantity C results from the evaluation of the trace over 
the internal atomic variables performed in the decoupled 
basis, where the dipole operator T> has simple matrix 
elements on 1/2 — > 1/2 and 1/2 — > 3/2 transitions. It 
only depends on the atom-laser detunings Sj: 



3 l + 25a/5i 



2 2 + Sl/Si 



(31) 



In principle, C can range between —3/4 and 3/2. In the 
useful case of detunings much larger than the fine struc- 



one has S 



3/2 



§1/2 and C approaches 



ture w 3 / 2 
3/2. 

The expression ([30]) allows a straightforward numerical 
evaluation of the heating rate. It is convenient to intro- 
duce the lattice depths Uq^ > along each direction 
fi G {x, y, z}, such that Eq. (fTB| reads 

C/ rod (r) = -P^C/o^cos 2 k L r„ (32) 
C4iuc(r) = P g ^2 U ^ sin 2 k L r^ (33) 



for a red detuned and a blue detuned lattice, respectively. 
In Fig |3l we plot the numerical results for the ground 
band heating rate in an isotropic lattice Uo ifi = Uq V/x, 
either red or blue detuned, in units of the total heat- 
ing rate (|25p . as a function of the lattice depth in units 
of the recoil energy ([9]), taking C = 3/2. Whereas the 
ground band heating rate is of the same order as the total 
heating rate for shallow lattices Uq < E tec , it drops to 
much smaller values in the tight-binding limit Uq 3> P re o 
the main effect being that the ground band width drops 
exponentially fast with (Uq/ E^) 1 / 2 in that limit. Fur- 
thermore, it is observed in Figj3]that the ground band 
heating rate drops faster for a blue detuned lattice than 
for a red detuned one. This we now explain analytically, 
not restricting to an isotropic lattice. 




FIG. 3: Initial ground band heating rate dEGs/dt for an atom 
prepared in the ground Bloch state of the optical lattice, for 
a red lattice (red lines) and for a blue lattice (blue lines), as a 
function of the lattice depth in recoil units, Uo/E Tec . Here the 
coefficient C = 3/2, and the lattice is isotropic: Uo tli = Uo for 
all spatial directions n = x, y, z in Eq. (|32[) for a red lattice and 
in Eq. (|33[) for a blue lattice. Solid lines: Numerical evaluation 
of (I30[) . Dashed lines: Asymptotic expressions ()40|) and (|42|) . 
The ground band heating rate is expressed in units of the 
total heating rate (|25[) . 



Lamb-Dickc regime: We first rewrite the matrix element 
appearing in pop in terms of the ground band Wannier 
function $o(r), assumed as usual to be real and normal- 
ized to unity: 



(0;fc L (e p - n)|e 



— ik,L n-r 



^W|0;0) = 



^2 [ d 3 r $ (r - R) e -^ n ' r £ M (r)$o(r) (34) 

R ^ 

where the sum is taken over the locations of the lattice 
potential minima R 6 (jr/k^Z 3 . In the large lattice 
depth limit Uo !f j./E ICC ^> 1, the ground band dispersion 
relation may be approximated by the tight-binding ex- 
pression 



Eo-a — Eo.O 



2t v [l - cos(q„ir/k L )}. (35) 



The sum (|34[) . to zeroth order in the tunneling amplitudes 
t v > 0, reduces to the R = term. This term may be 
evaluated by a Lamb-Dickc expansion in powers of /clo m , 
resulting from series expansions in powers of k^r, e.g. 



1 - ik L n ■ r + 0{k L rf 



We have introduced the size 

1/2 



2niui, 



1 / E re 



k L \ 4 ^o,p 



1/4 



(36) 



(37) 



of the ground state of the harmonic oscillator of angular 
oscillation frequency u>^ that approximates U(r) around 



r = along direction \x. Restricting to leading order in 
the Lamb-Dicke expansion, we can also approximate $o 
with the Gaussian wavefunction of the harmonic oscilla- 
tor ground state. 



Atom 


Wavelength Xl 


E rec 






Ref. 


87 Rb 


850 nm 


3.2 kHz 


6 s" 1 


1 s- 1 


[22] 


87 Rb 


750 nm 


4.1 kHz 


12 s- 1 


2 x 10- 2 s- 1 




133 Cs 


1064 nm 


1.3 kHz 


1 s" 1 


0.1 s- 1 


[23] 


133 Cs 


780 nm 


2.5 kHz 


2 s- 1 


5 x 10 -3 s- 1 




87 Rb 


x axis: 765 nm 


3.9 kHz 


8 s- 1 


10- 2 s- 1 






y+z axis: 844 nm 


2.5 kHz 


2 s- 1 


0.3 s- 1 






Total: 




10 s- 1 


0.3 s- 1 


m 



TABLE I: Total and ground band heating rates for two dif- 
ferent experiments with 87 Rb and 133 Cs, for blue and red de- 
tuned isotropic optical lattices of depth per axis Uo,^ = 10 En. 
The heating rates are estimated using the formulas (|40|l and 
(|42p obtained in the Lamb-Dicke limit, and normalized by the 
ground band tunnel energy t„ ~ 2 x 10- 2 -Br defined in (|35|) . 
For these parameters, the coefficient C in Eq. (|3ip is within 
10 % of the asymptotic value 3/2. The corrections to the light 
shift potential U discussed in subsection III CI are very small, 



1Q- 



10" 



SU n - 



10" 



10-', 5U 1 



/e 2 



10" 



The last two lines refer to the bichromatic lattice of || 

Red lattice: We use the expansion £ M (r) = £o,n + 
0(/cl^j) 2 so that 

(0; k L {e^ - n)|e-^"-%(r)|0; 0) = £ Jl + 0(k L a ho ) 2 }, 

. (38) 

where a\ w is the largest of the a v . Upon insertion in 
pi!]), we face angular integrals that can be calculated in 
spherical coordinates of axis is, e.g. 

f d 2 n 2 f 1 du o , , 2 

I n cos(7T7i^) = / — u cos(7ru) = -. (39) 

J 47T J_ x 2 7H 

We obtain the deep lattice asymptotic expression at time 
t = 0: 



d 2 j 



fl,U 



C 

— (1 + V) 



(40) 



where 8^ is the Kronccker symbol. Expression (|40[) 
asymptotically matches the numerical results, see Fig|3l 
We thus find that, for a red detuned lattice, the ground 
band heating rate is smaller than the total heating rate 
by a factor scaling as the tunneling amplitude over the re- 
coil energy. This reduction factor simply originates from 
the energy width of the ground band. 
Blue lattice: The laser field amplitudes vanish in r = 
and have the expansion £^(r) = S^^k^r^ + 0(ki J r fJ ) 3 so 
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that [Hi 



(0;fc L (e M - n)|e 



^W|0;0) = 



£o, /4 H(fcia /1 ) 2 +0(fc i a llo ) 4 ]. (41) 

Insertion in p0[) and angular integration give at time 
t = 0: 



dt 



E 



blu< 
GB 



^2^|£ , Al | 2 (fcLa Al ) 4 



3 7T 

i 

— + 



33 

^4 



5 7T 
11 

3^2 



(42) 



This can be further simplified, since |f o,/j, | 2a f t is actually 
independent on the amplitude of the standing wave 
and hence on the direction of space. Expression (|4"2")l 
asymptotically matches the numerical results, see Fig[3l 
We thus find that the blue detuned lattice has a ground 
band heating rate that is reduced as compared to the 
red one (|4"2"j) by factors (/cl^i) 4 oc E lec /Ua !fl <C 1. The 
first factor (fc^a^) 2 originates from the reduction of the 
fluorescence rate as compared to the red lattice, due to 
the Lamb-Dickc effect: Absorption of a laser photon in 
the standing wave £o,fi sin fe^ry, brings the atom mainly in 
the first excited band, with a small transition amplitude 
oc kLdft leading to a reduced fluorescence rate oc (fc^a^) 2 . 
The second factor (^lo^) 2 originates from the branching 
ratio of spontaneous emission from the internal state e 
in the first excited band to the internal state g in the 
ground band. 

Limits of validity: The ground band energy increases lin- 
early in time (with the rate that we have calculated) 
only for short times. We can estimate the onset of 
non-linearity by considering double fluorescence cycles 
(with the atom arriving in the ground band after the 
second cycle). For a red detuned (respectively blue 
detuned) lattice, the probability of such a double cy- 
cle over the time interval t is cx (rg^i) 2 [respectively 
oc (rg^£) 2 (fc£,aho) 4 ], corresponding to a ground band 
energy increase oc t v . The corresponding mean increase 
of energy is negligible as compared to the single cycle con- 
tribution r^tt v [respectively T^tt u (k L a ho ) 4 ] as long 
as r™ u a ^t < 1 in both cases. 

As a final remark, we point out that the ground band 
heating rate for a blue lattice, being much smaller in the 
tight binding limit than the one for a red lattice with the 
same depth Uq/ 'E Tec , is thus much more sensitive to any 
additional heating effects not included in our theoretical 
model, in particular to collisions between ground band 
and excited band atoms in the many-body case. 



V. CONCLUSION 

We have performed a full quantum calculation of single 
atom heating rates in a far-detuned, three-dimensional 



optical lattice, including explicitly a realistic atomic in- 
ternal structure and the fact that the superimposed laser 
standing waves have different frequencies, as in real ex- 
periments. 

First, we have calculated the total heating rate, that is 
the rate of increase of the total mechanical energy of the 
atom. Remarkably, we have found a universal expression, 
independent of the initial internal and external atomic 
state, and simply equal to the product of the recoil en- 
ergy E rec and of the maximal fluorescence rate Tg 1 ^ that 
may be realized in the lattice. The total heating rate is 
thus independent of the sign of the laser detuning (red 
or blue). 

This general feature is easy to understand in the lim- 
iting case of a deep lattice for an atom initially in the 
ground energy band: The total heating rate is then de- 
termined by rare photon scattering events transferring 
the atom to excited bands at a rate which is indepen- 
dent of the sign of the detuning. In the blue-detuned 
case, the atom most probably arrives in the first excited 
band after a scattering event, which increases the en- 
ergy by the oscillation quantum Huj OSC7 that is by much 
more than the recoil energy; this however takes place at 
a rate oc r^ u a *i? rcc /fiw osc <C r^f* because of the Lamb- 
Dicke effect. The product of the rate and of the energy 
change is indeed oc Tf^*E rcc . In the red-detuned case 
[I3|, the fluorescence rate, of order Tg 1 ^ , is much larger; 
the atom however mainly returns to the ground band 
after a scattering event, due to the Lamb-Dickc effect, 
which increases the energy by much less than the recoil 
energy; with a branching ratio oc E Tec /hio osc , the atom 
however arrives in the first excited band, which increases 
the energy by hu osc and results in a heating rate again 

^ 1 fluo "^rec- 

Second, for an atom initially prepared in the lowest 
Bloch state of the lattice (a minimal single-particle model 
for the many-body superfluid state realized in experi- 
ments), we have calculated the initial ground band heat- 
ing rate. This is the rate of energy increase due to pro- 
cesses where the atom returns to the ground band of the 
lattice after a photon scattering event. We have derived 
analytical expressions of this rate in the deep lattice limit. 
They show that, in contrast to the total heating rate, the 
ground band heating rate strongly depends on the laser 
detuning, and is strongly suppressed for blue detuned lat- 
tices: It is of order rg^to f° r a red deep lattice, and of 
order T^^to(E ICC /huj osc ) 2 for blue deep lattice, where to 
is the atomic tunneling amplitude between neighbouring 
sites. 

A recent experiment Q reported a measured heating 
rate of (dE/dt) oxp /to ~ 3.5s -1 , where the tunnelling am- 
plitude to was essentially independent of the spatial di- 
rection. This heating rate was obtained by using Quan- 
tum Monte Carlo simulations at several temperatures to 
calibrate the experimental data: The decrease of the 
visibility of the interference pattern observed after free 
expansion was recorded over time, and compared to a 
"visibility-energy" abacus obtained from the simulations. 
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The heating rate measured this way is unsurprisingly 
smaller than the total heating rate that we calculate, 
see Table HI Indeed, the Quantum Monte Carlo simula- 
tions take into account the ground band only, and we ex- 
pect the measured interference pattern to depend mostly 
on the ground band atoms. Hence, the most pertinent 
rate to compare the experiment to is the ground band 
heating rate. The measured heating rate is significantly 
larger than the ground band heating rate, though, see 
Table HI To resolve this discrepancy, one would have to 
include effects that we did not consider in this article: 
Firstly, heating due to technical noise in the apparatus, 
which should be quantified for a specific experiment, and 
secondly, the role of collisions in redistributing part of 
the energy from the excited to the ground band. To 
our knowledge, the latter, more fundamental many-body 
problem, is still quite unexplored [24j and provides an 
interesting direction for future work. 

We acknowledge useful discussions with I. Bloch, S. 
Trotzki, L. Pollet, N. Prokof'ev, B. Svistunov, M. Troyer, 
K. Mur, and J. Dalibard. The authors are members of 
IFRAF. This work was supported by the DARPA project 
OLE. 



Appendix A: Derivation of the ground state master 
equation 

In this Appendix, we provide details about the deriva- 
tion of the master equation used in the main text. We 
use an interaction picture with respect to the kinetic plus 

2 

hyperfine Hamiltonian Hq = j*— + H^f, where operators 

arc modified as X(t) = e iHot / h X(t)e- iHot / h . The Bloch 
equations read 

r 



-a ejg = (iSj - -)a ej 

e.ig 



ill 



V / era \ 



e j'9 



(Al) 



state optical coherence 



cr 



"tL^-myr (v AL (t-r) 



ih 



a gg (t-T). 
(A4) 

By repeated integrations by parts of the integral over 
t, integrating the exponential factor, we get a formal 
expansion of a ejg (t) in powers of 1/(5 j + iT/2), that we 
turn into a series in 1/5 j. 

To order 1/5: The optical coherences in Schrodinger pic- 
ture are given by 



1 



= Wj (VAL(t)) ejg a gg (t). 



(A5) 



This is familiar for a constant atom-laser coupling. We 
have shown that it holds even for a time dependent cou- 
pling provided that Eq. ([3]) holds. Inserting this relation 

Eq. (|A5|) in Eq. (|A2| gives in steady state a^} ejl = 0. 
After insertion of Eq. (|A5|) in Eq. (|A3|) we then find that 
<7gg ag has a purely conservative evolution of Hamiltonian 

ffjjf + + U^(r, t) with the time-dependent light shift 
potential: 

=E^j( y AL(t)) flej (V AL (t)) ejg . (A6) 

As expected U^/h is of order VI 2 / 5 so that dd gg /dt is at 
most of this order; in general, can induce Zeeman 
and even hyperfine ground state coherences (32| 
To order 1 /5 2 : From Eq. (|A4[) we obtain 

hM^) e jM- (A7) 



We shall neglect da gg /dt with respect to Ta gg as allowed 
by Eq. ©. Inserting Eq. (TAT]) in Eq. ([A"3j) has two effects. 
First it induces a small modification of the light shift 
potential of order Eq. ([H|) : 



U 



(2) 



f 



X/ (t 



(fr5j) 
ih[ , rr . „„ ( d 



j[(V Ah (t)) ge (-V Ah (t) 



-h.i 



}. (A8) 



dt aeje " 
1 

ih 



dt 



V A h{t)\ a gej , - a ejg (ValW 



(A2) 



ifi? ( 



Val(*) 

ge.J 



h.c. 



(A3) 



The feeding term of the ground state by spontaneous 
emission is given in Schrodinger picture by ([TB")) . 

We perform the series of approximations discussed in 
the main text. Integrating formally Eq. (|AI[) . after omis- 
sion of (J ee , neglecting p 2 /2m as compared to HT, and ne- 
glecting a transient of duration 1/r, we obtain the steady 



Second it induces a lossy evolution of cr^ 8 , Jjc g£ , 



-\{W,(X gg }, with 

W (2) = y 1 



d 

T+ dt 



(VAL(t)) gej (V Ah (t)) ej 



9 ' 



(A9) 

The adiabatic elimination of the excited state matrix el- 
ements up to order 1/5 2 , remarkably leads to an expres- 
sion similar to the usual two-level atom case, despite the 
presence of the hyperfine Hamiltonian and the time 
dependence of Val: 

(Val(*)) 



diag 



7 (2) 
'ejeji 



ejg a gg 



!l< J 



h 2 5.j5 



-+0 



joj> 



9l 

5T 
(A10) 
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This may me checked by direct insertion of Eq. (|A7[) in 
Eq. (|A2|) . the 0[. . .} term coming from neglecting da gg /dt 
as allowed by Eq. (|8|). When inserted in Eq. (|A3|> this 
provides a closed equation for cr^ ag . 
Temporal average: Under the conditions discussed in the 
main text, we now average out the rapidly oscillating 
terms in the ground state master equation. The lattice 
potential U^'(t, t) to order 1/8 contains a zero frequency 
part, which is scalar according to Eq. §5§ and is called 
U in Eq. (|T5|) . It also contains oscillating contributions 

C/^ ) (r)e" i;A ^ t , with A^ = - A„ ^ 0, fx ? v : sup- 
posed to be pairwise distinct, as it is the case in and 



(i^) gej £*(r)(Z^) ejg £ M (r) 

Mr 



(All) 



These contributions are in general not scalar. It is con- 
venient to split l/jh) in (i) an off-diagonal part, inducing 
coherences between the two hypcrfine ground states, that 
is a raising (resp. lowering) part (ufu))± coupling cji to 
g 2 (resp. g 2 to 51), and (ii) a diagonal part (£/$/)diag 
coupling gi to g\ and (72 to g 2 . In the interaction pic- 
ture, the diagonal parts correspond to terms modulated 
at frequencies A^ in U^-\ and the off-diagonal parts to 
terms modulated at frequencies A M „ =p Aj^ f . They will 
induce rapidly oscillating components of the density op- 
erator Ggg at those frequencies, these components being 
very small here since Ml 2 /8 <C lA^j, A^ f . In the present 
regime Af lf 3> A mo d, we can treat separately the effect 
of the off-diagonal and diagonal parts of U^ 1 '. 

The off-diagonal part of U^(t) induces a Rabi cou- 
pling ~ HQ, 2 1 '5 between g\ and g 2 much smaller than 
their energy splitting hAf l{ , with a very slow time varia- 
tion at the scale of 1/A^ f . Hence we treat it directly to 
second order in usual perturbation theory for a fixed time 
t and then average the result over time. This leads to an 
effective light shift potential acting within each hypcrfine 
ground state: 



(i) 



off, I 



J2 (uftUufth (uft)+(u$)- 



(A12) 

This corresponds to Eq. ([TT]) with a different notation. 
Having eliminated the hyperfinc coherences, we consider 
in what follows o-g g as - 

The diagonal part of Uj^J couples the zero fre- 
quency component 5^ ag to rapidly modulated com- 
ponents <j£g (t)e~ tA ^ t ; adiabatic elimination gives the 
slowly evolving amplitudes 



[(^)diag,^] 



hA 



(A13) 



//// 



They are small according to Eq. (fT0| . Including the cou- 
pling of af g as to o££ by (UfiJ) dias gives the a priori lead- 
ing contribution to the effective potential induced by the 



diagonal micromotion: 



U. 



(i) _ NT 

off, II ~ / j 



[(^M^)diag' (^w)diag] 



2HA, 



(A14) 



Usually, the micromotion is studied for a spinless parti- 
cle, in which case the commutator in Eq. (|A14[) automat- 
ically vanishes and the first non-zero correction scales as 
l/A^ od [13, [ll|. Here, the atom has a non-zero s pin 
From a calculation in the decoupled basis we find [25( 



U, 



(i) 



{2id 2 /3h 2 ) J ■ (e M x g v )[8-} 2 -5-}± The com- 



^ = £ ll £*(2id 2 /3ti 2 )J-{e f _ l xe„)[S 3/ 1 2 -v 1/2l 
mutator in Eq. (|A14|) thus vanishes also in our spinorial 
case for the field Eq. (HI). G oing to next order in 1/ A mo d 
with the formalism of [2fJ extended to the case with a 
spin we finally find a non-zero contribution 



V. 



(i) 



grad(C/W) dia 



?rad (Z7$/)dia« 



cff.II ~~ / j 



2mAJ„ 



(A15) 



This corresponds to Eq. (fl3|) with a different notation. In 
appendix [B] we present an alternative derivation of this 
result, based on the dressed atom approach, and showing 
that there is no other micromotion term of order 1 / Aj^ od . 
Summing both contributions Eqs. (|Al2|Al5j) . we get the 
effective time-independent correction to the light shift 
potential for the theory of order 1/6, integrating out the 
effect of hypcrfine couplings and diagonal atomic micro- 
motions. 

To eliminate the rapidly varying temporal components 
in the 1/8 2 terms of du gg /dt, we simply take the temporal 
averages, neglecting micromotion corrections at this or- 
der, as allowed in particular by Eq. ©. We also project 
out the operators inducing ground state hypcrfine co- 
herences. The 1/8 2 correction E^jfag t° the light shift is 
obtained by projecting out the hyperfine coherences of 
Eq. (|A8|) and taking the temporal average. Its explicit 
expression is not required here. The complete expression 
for SU in Eq. (fT5]) is 



SU = u 



ff.I 



u 



(1) 



off, II 



u 



(2) 
diag' 



(A16) 



The time average of Eq. (|A9[) and the use of the scalarity 
relation Eq. ([3]) gives the lossy part of Eq. (|15p . The 
time average of Eq. (|A10[) gives Eq. (|19j) with a different 
notation. 



Appendix B: A derivation of the micromotion 
potential based on the dressed atom picture 

We propose here a derivation of the micromotion po- 
tential alternative to [20j and including the atomi c sp in. 
The idea is to use the dressed atom approach |26| to 
eliminate the time-dependence of the Hamiltonian, and 
to use standard time-independent effective Hamiltonian 
theory [2(|. The laser field is then assumed to be ini- 
tially in a Fock state with huge occupation numbers 
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in each mode of frequency ojl + A M . One may then 
neglect the dependence of the atom-laser coupling with 
the photon number, replacing there the photon annihila- 

1 /2 

tion operators a M with A^, where the phase operator 



basis 
with 



[a u a 



V) 1 ^ 2 a^ has unit matrix elements in the Fock 
The resulting Hamiltonian is H = Hq + V 



H = fto + X) fi ( Wi + A /*)(°l» f -n A*) ( B1 ) 



(B2) 



where /i = E— + U is a purely atomic Hamiltonian, 

C/^y is given by (|A11|) and ?7 is the time independent 
part of U {1 \ called U in ((161). It is convenient to set 

(f^i )diag = Zt^pTfu, where the operator does not 
depend on the atomic position and obeys 



rp\ _ rp _ rp 



Since the field is initially in the Fock state 
\( n n)nE{x,y,z})i we introduce the orthogonal projector P 
on that Fock state, and Q = 1 — P the supplementary 
projector. The effective Hamiltonian H c s(E) inside the 
subspace over which P projects is exactly given by 



H eS (E) = h P + PVQ 



Q 



EQ - QHQ 



QVP, 



(B4) 



where we used PVP = resulting from the fact that 
V does not conserve the number of photons in each 
mode. At this stage, E is arbitrary but much smaller 
than ftA mo( j. Since the denominator in (|B4[) is of order 
fiA mo( j » V, because of the occurrence of the projector 
Q, we may expand (|B4[) in powers of V, restricting to 
terms of order up to l/A^ od : 



H eS (E) = h P + PVQ 



PVQ 



EQ - QH Q Q 

Q Q 



QVP 



EQ - QH Q Q 



QVQ 



EQ - QH Q 



QVP 



O(l/A 3 mod ). (B5) 



We first focus on the terms of order V 3 in (|B5[) . In 
the denominators, H(uil + A^)(aj l a fl — n M ) gives con- 
tributions of order fiA mo d that dominate over E — ho- 
Expanding in powers of (E — ho)/A Tao d then gives 

H eS (E)- ^ P\£ a \ \£p\ \£,\ [ h 2 AagA _ 



h 2 A a pA ai 
where we recall that A„ v 



P + 0(l/A 3 mod ), (B6) 



A v , and the prime 



pairwisc distinct. Since the modulation frequencies A^ 
are incommensurate, a first action of V on the Fock state 
|(?V)), e.g. creating a photon in mode a and annihilating 
a photon in mode /3, has to be compensated in two steps, 
either annihilating a in the first step and creating (3 in 
the second step, or vice-versa, hence the occurrence of 
two terms in (|B6|) . Exchanging the summation indices 
a and (3 in the second term of (|B6|) . and using the anti- 
symmetry of T a p and A Q( g under the exchange a <-> /?, a 
sum 1/A^ 7 + 1/A 7| g = appears, so that 



H^(E) = 0(l/A 3 ulod ), 

and may be neglected. 

We are left with the V 2 contribution to H c ff. 



(B7) 



(B8) 



(B3) We expand this expression up to order 1/A 



mod ' 



H$(E) = ^P|^| 2 |^| 2 T M ,— \-T vll P 

+ O(l/A? nod ). (B9) 



{hA^Y 



The first contribution, of order l/A mo d, corresponds to 
the effective potential given in (|A14[) and vanishes, as al- 
ready noted in Appendix ^ since and A M „ are anti- 
symmetric functions of y, and v. The second contribution 
does not vanish. Since it is of order 1 / A^ od already, we 
treat it in perturbation theory: An unperturbed eigen- 
state | Ao) of ho of energy Eq experiences an energy shift, 
calculated here up to order (]//A mo d) 2 , equal to 



E2 = (\o\H$(E Q )\\o) 



(BIO) 



Taking advantage of the fact that A^„ = A 2 , we use 
the relation 



(Ao|(^ ) )diag(/ i Q - £o)(C/W)diag|A ) + M <"> V = 

(Ao|[[(C/W)diag,M, (^^diagHAo). (Bll) 

Since U is scalar, it commutes with (£/^)di ag and the 
kinetic energy operator p 2 /(2m) is the only one to con- 
tribute to the commutator. After an explicit calculation 
of the double commutator, we obtain 



E 2 = <Ao|[/« n (r)|A 



(B12) 



on ^2 means that the sum is taken over indices that are where U^g U (r) is indeed given by (|A15[) . 
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